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ABSTRACT 

We present new V- and i?-passband CCD photometry of UU Sge and 
V477 Lyr, the echpsing binary nuclei of the planetary nebulae Abell 63 and 
Abell 46, respectively. We have performed a simultaneous analysis of VR light- 
curves and estimated the effective temperatures for the primary and secondary 
stars to be 78 000 ± 3000 and 6136 ± 240 K for UU Sge, 49 500 ± 4500 and 
3874 ± 350 K for V477 Lyr. We have also reanalysed the previously mea- 
sured radial velocities and combined the results with those obtained from the 
analysis of the light curves to derive absolute parameters of the components. 
The secondary stars have larger radii than expected from their main-sequence 
counterparts at the same masses. We have determined the post-common en- 
velope ages and the thermal time scales of the systems and examined the 
possible reasons of expanded radius of the secondary components, together 
with some selected post-common envelope binaries. We conclude that the 
secondary components of the nuclei of the planetary nebulae are still out 
of thermal equilibrium along with two post-common envelope systems: HS 
1136+6646 and RE 1016-053. For other systems, magnetic activity has been 
suggested as the more plausible reason for their expanded radii. We have also 
estimated the common-envelope efficiency parameters of UU Sge and V477 
Lyr. 
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1 INTRODUCTION 

The classical approach to planetary nebulae (PNe) is that they are envelopes ejected by 
single stars (0.8^M/Mq^10) at the end of the asymptotic giant branch (AGB) stage. How- 
ever, up to date, it has been reported that at least 10-15 per cent of planetary nebula 
nuclei (PNNi) are very close binaries with periods in the range of hours to a few days 



(IBond fc Grauer 



1987 



Bond &: Livio 



199' 



Bond 



2000) • Recently, some radial velocity sur- 



veys have been carried out in order to searc h for new variab 



long to be detected using photometry (e.g. 



De Marco et al. 



e FN 



2004 



Ni, which have perio ds too 



Afsar fc Bond 



20051). The 



existence of these o bjects has 



by several authors (iPaczyhski 



)een e xplained by means of common enve. 



1976 



Iben fc Livio 



1993 



Iben fc Tutukov 



ope (C E) ejection 



19931 ). According 



to CE evolution theories, the atmosphere of the primary component will expand during the 
red giant branch (RGB) or AGB evolution, and depending on the initial separation of the 
binary system, both components will be engulfed by extended atmosphere of the primary. 
The gravitational drag forces arise due to orbital motion of the secondary component around 
the primary's core (spiral-in process) will cause orbit to shrink dramatically. Because of the 
conservation of total angular momentum, the orbital angular momentum will be transferred 
to the CE and as a result the envelope is spun up and ejected from the system. Thus, the 
outcome will be a close binary as a nucleus of a planetary nebula. The efficiency of which 
orbital energy is deposited into CE ejection is defined as "efficiency parameter", denoted as 
acE- Our knowledge on the CE evolution is highly limited due to the difficulties in deter- 
mining of stellar and nebular parameters. However, the central stars of the PNe, which are 
also members of eclipsing binary systems, are the direct evidence of the existence of the CE 
phase and provide an opportunity to drive accurate absolute parameters of these objects. 
There are a number of these objects and only four of them are eclipsing binaries: UU Sge, 



V477 Lyr, BE UMa and SuWt 2. ' 
either due to reflection effect (see 



'he other close 



)inary nuclei are known to be as binaries 



De Marco et al. 



20081 for a detailed study) or ellipsoidal 



variation they show in their light curves, or they are just single-lined spectroscopic binaries. 

The eclipsing binary nuclei UU Sge, the nucleus of Abell 63, and V477 Lyr, the nucleus 
of Abell 46, were previously studied by several authors. The eclipsing nature of UU Sge 



was flrst reported by Hoffieit in 19 32. A 



brightness nebula was discovered by 



t er th ree decades the pr esence of a low-surface 



Abelll (Il966l ). and iBondl (119761 ) showed the connection 
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between PN and UU Sge . The first photoelectric photometric observations of the system 



were made by 



Bond et aL 



(119781 ). They obtained a hght curve in B- passband and indicated 
that the system has a total eclipse at the primary conjunction. They also noted the extreme 
heating effect in the system due to large temperature difference between the components. 
From the solution of the light curve they concluded that the system contains an sdO primary 
of ~0.9 Mq with a temperature of about 35 000 K and a dK secondary of ~0.7 Mq. They 
also emphasized that the observations of the system was contaminated b y the nearby optical 
companion ~4" (later to be measured as ~2".8 by IWalton et al.lll993l ) to the east of UU 
Sge and they could not determine the true depth of primary minimum accurately. However, 
they estimated the contamination of the optical companion from the primary minimum 
observations and suggested that the actual depth of the primary eclipse should be about 4.3 
mag. 
Later, 



Budding fc Kopall (jl980l ) obtained the F-passband light curve of the system. They 
found the 0-type primary component has a tem perature of not rnuch hi gher than 30 000 K 
and the secondary component is around 6000 K. IVan Buren et al.l (jl980l ) discovered the soft 
X-ray emission feature of UU Sge and suggested that the observed X-ray originates in the 
cor ona of the late-type st ar . 



Walton et al. 



(Il993l ) made spectroscopic observations of the system at several phases. 
Using the spectrum taken during the primary minimum they estimated the spectral type of 
the secondary component as G7 V. From the spectra obtained at out-of-eclipse they deduced 
that the exciting hot central star is an sdO type PN nucleus with a temperature of around 
50 000 K. 



The first V- and /-passband CCD photometry was made by lPoUacco fc Belli (11993. hereafter PB93I ). 
However, in /-passband, they only obtained the primary minimum since they aimed to use 
the primary eclipse to determine an accurate temperature for the secondary. Using the 
medium-resolution spectroscopy, they obtained the first radial velocities for the system. 
They used the results deduced from the analysis of radial velocities for the solution of the 
F-passband light curve. Since CCD observations enabled to subtract the contamination of 
the neighboring field star, they derived the more accurate absolute parameters for the in- 
dividual components. After performing a series of light curve solution, they estimated an 
effective temperature of 117 500 K for the primary component and fixing this value for 
the final light curve solution they suggested a temperature of ~7300 K for the secondary 
component. 
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In a later study, iBell et al. 



(11994. hereafter BPH94I ) confirmed that the primary eclipse 
is total with a duration of ~14 min. They also obtained the spectroscopic observations of 
the syste na during primary minima and extracted the spectrum of secondary component. 
Applying iKuruca (Il979l ) atmosphere models to the spectra they determined a secondary 
component's temperature of 6250±250 K. Adopting this value for the secondary component 
temperature they revised the preliminary light curve analysis. Solutions ind icated a primary 



component with a temperature of 87 000±13 000 K for a limb-brightened (iBond &: Grauer 



19871 ) secondary component. 



V477 Lyr, the eclipsing binary nucleus of the PN Ab ell 46, has a similar light curve 
to UU Sge, except that the eclipse is partial ( Bond! 1980). The system shows a very large 
reflection effect which creates sine-like distortion on the light curve with an amplitude of 
~0.6 mag. 

Although its eclipsing nature was reporte d more than two decades ago, a few light curve 
solutions of V477 Lyr have been published. iRitten (119861 ) presented a solution based on a 
light curve analysis by L.W. Twigg (communicated by Bond, 1981). This solution yielded 
an orbital inclination of 80?8 with a mass ratio of Mi/M2~2.3, radii of -Ri~0.09 and 
i?2~0.20 Kq. The temperature of the primary component was taken as Ti ~105 000 K 



Cons idering a typical mass value for the central star of a planetary nebula ( ISchonberner 



19811 ). Ritter assumed a mass for the primary cor aponent as Mi~0.6 M p), which yields a 



secondary component mass of M2~0.25 M0. Also, iBond &: Grauerl (119871 ) cited the results 
of a preliminary analysis by J. Kaluzny and a temperature of 60 000 K for the primary 
component was suggested. 

The most recent and detailed solutions of both light and radial velocity curves were made 



by iPoUacco fc Belli ( 1l994. hereafter PB94I ). They obtained the medium- resolution spectra 



and a new F-passband CCD photometry of the system and presented the first complete 
light curve of V477 Lyr. They performed a series of light curve solution for the effective 
temperature of primary component ranging from 30 000 to 130 000 K. Using the mass 
ratio (M2/Mi=0.29) obtained from the radial velocities and adopting the limb-brightening 
for the secondary component, they achieved a best fit to the observations for the effective 
temperature of 60 000 K. Thus, the F-passband light curve solution yielded an inclination 
of ~80?5, and a temperature of ~5300 K for the secondary component. 

In this paper, we present new two-colour CCD photometry of UU Sge and V477 Lyr. 
The aim is to constrain the effective temperatures of the components through simultaneous 
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solutions of the two-colour light curves. We will estimate the common envelope efficiencies, 
acE, and discuss reasons for the oversized secondaries that we find. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 UU Sge 

We have obtained the Johnson V- and i?-passband photometry of the system at the National 
Observatory of the Scientific and Technological Research Council of Turkey (TUBITAK- 
TUG) in Antalya. TUG is established on the mountain of Bakirlitepe (in the West Taurus 
Mountain ranges), at a height of about 2500 m. Observations were carried out using a 1.5-m 
Cassegrain telescope (RTT 150) equipped with 2048x2048 back illuminated CCD camera 
which has an area of 9'.lx9'.l and a scale of 0.26"/pixel. We observed the system on the 
nights between July 4 and 7 and on the nights of 2002 August 6, 7 and 11. Since the orbital 
period of UU Sge is relatively short (~11.5 hours) and the system has a very deep primary 
minimum (~4.5 mag), exposure times were li mited between 60 and 120 s. The light curves 



(120001 ) 



were phased by using the ephemeris given by iKiss et al 
HJD (Min.I) = 2451766.5285 + 0.46506921E d. 

UU Sge has a nearby optical companion ~2".8 to t he east of the b i nary nucleus. The 



visual magnitude of the companion is V=15.87 mag. (ICiarduUo et al.l Il999l ). The visual 
magnitudes of the binary system itself at maxima and in mid-eclipse are Vmax=14.67 mag 
and Vrrtm=19.24 mag (PB93), respectively. 

The data were reduced using the appropriate tasks in irafQ. First a master bias was 
created and subtracted from the images. Then we used an averaged master flat-field, which 
we created using the flat frames taken both during the evening and morning twilight, and 
normalized all the images. Since UU Sge is located within a crowded region, images were also 
trimmed right after bias and flat-field corrections. After ffist reduction steps were completed, 
iRAF/DAOPHOT package was used to apply the crowded field photometry. Thereby, we were able 
to make profile fitting and PSF (Point Spread Function) modeling of the stars in the target 
field and subtract the light contamination of the nearby optical companion. We measured 
the PSF magnitudes of UU Sge and five nearby stars to apply differential photometry. These 



^ IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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Figure 1. V- and fl-passband light curves of UU Sge (filled and open circles, respectively). 

five stars were chosen as comparison star candidates which were previously determined by 
PB93 (CI, C2, C3, C5 and C6: Fig.l, therein). After investigating all comparison candidates 
for a possible variability, wc decided to use C3 as the comparison, and we determined the 
differential magnitudes for UU Sge as plotted in Fig.l. 

We have determined the depth of the total eclipse to be 4.450±0.075 and 3.764±0.074 
mag for V- and i2-passband, respectively. The depth for the F-passband is somewhat 
larger than that of BPH94. The time difference between the fourth and first-contact points 
was determined as ~1.26 hour, and the duration of the total echpse at the bottom of the 
primary minimum was estimated as ~13.4 min. Not only the short time length but also the 
integration time increase at the bottom of the total eclipse, with a relatively long uploading 
time of the camera, brought about to obtain only a few data points in conjunction. Moreover, 
the minimum brightness of UU Sge is close to the magnitude limit (in case of photometry) 
for a 1.5-m telescope, which led to have scattered data points at the primary conjunction 
with errors around as much as 0.06 mag. 

2.2 V477 Lyr 

Two-passband (Johnson V and R) CCD photometric observations of V477 Lyr were carried 
out using the same setup as described for UU Sge. We observed the system during the nights 
from 2002 August 6 to 10. V477 Lyr has a partial eclipse at the primary minimum and it is 
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Figure 2. V- and ij-passband light curves of V477 (filled and open circles, respectively). 

not as deep as that of UU Sge. We set the exposure ti me as 60 s during the overall phases. 



Kiss et al 



feood), 



The ephemeris used to phase the data was taken from 

HJD (Min.I) = 2451766.4038 + 0.47172909£' d. 

After bias subtraction and fiat-fielding, the data reduction was performed applying the 
aperture photometry package of iraf/apphot. The magnitudes of the stars including V477 
Lyr and other five comparison candidates, which were previously found to be constant by 
PB94 (CI, C3, C4, C5 and C6: Fig.l, therein), were measured. The errors of the magnitude 
measurements at primary minima did not exceed 0.018 and 0.010 mag in V- and i?-passband, 
respectively. The differential magnitudes were computed with respect to the comparison star 
C6. The light curves of the system obtained in V- and /Z-passband are shown in Fig. 2. 



3 REANALYSIS OF THE RADIAL VELOCITIES 

The first radial velocities of UU Sge and V477 Lyr were previously obtained and published 
by PB93 and PB94, respectively. By using both absorption and emission lines (which move 
in anti-phase with the absorption lines and thought to originate on or very close to the 
surface of the secondary component) they measured the radial velocities of the systems. 
They used one of the spectra of the related object (UU Sge or V477 Lyr) as a template to 
calculate the cross-correlations. Therefore, as mentioned by the authors, only the estimates 
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Table 1. Results of the radial velocity analysis. 



Targets 


Ki 


K2 


q 


asini 




(km,/ s) 


{km/ s) 




{«0) 


UU Sge 


83.8 


182.8 


0.458 


2.45 




±6.7 


±7.4 


±0.034 


±0.09 


V477 Lyr 


51.9 


182.1 


0.285 


2.18 




±6.4 


±7.1 


±0.026 


±0.09 



of the velocity semi-amplitudes are meaningful while the systematic velocities of the systems 
remain unknown. That is why we were not be able to make simultaneous light and velocity 
curve solutions. However, we reanalysed the radial velocities by using the corrected radial 
velocity data given by PB93 and PB94 in order to increase the precision of the parameters 
and used them as inputs during the simultaneous solution of the two-colour light curves. The 
revised results of the radial velocity curve solutions are given in Table 1. The parameters Ki 
and K2 stand for the velocity semi-amplitudes of the primary and secondary components, 
q is the mass ratio (M2/M1) and a is the projected semimajor axis of the related system's 
orbit. 



4 SIMULTANEOUS LIGHT CURVE SOLUTIONS 
4.1 UU Sge 

We have performed a simultaneous solution of the V- and -R-passband lig 



empl oying the updated (2003) version of Wilson-Devinney (jWilson fc Devinney 



it curves by 



1971 



Wilson 



I979I . 1990, 2003) code (hereafter W-D). This program is written for the analysis of light and 



velocity curves of binary stars which have either circular or eccentric orbits. The program has 
different 'modes' in terms of Roche lobe filling conditions of the components. Examination 
of the light curve of UU Sge revealed a detached system and a circular orbit as expected for 
such a close binary. Therefore, we used 'Mode 2' which is developed for detached binaries. 
Since the reflection effect is highly strong as seen in the light curves, a detailed treatment 
to the reflection effect was applied during the solutions. 

The close binary systems, such as PNNi or post-CE (post-common envelope) binaries 
(subdwarf or white dwarf (WD) primary + low mass secondary) have large effective temper- 
ature differences between the components and show extreme reflection effects in their light 
curves due to irradiation around the substellar points of the cool secondaries. The modeling 
of the light curves of such systems is quite difficult and requires special treatment. As dis- 
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cussed by 



Hilditch et al. 



(119961 ). the temperature on the irradiated hemisphere of the cool 
secondary star is highly enhanced, therefore, contrary to the expected values (^42=0. 5, as 
generally assumed for cool stars) an albedo around or greater than one ( A2^^) is needed fo r 



securing a good fit to the detailed shape of the reflection effect (see also lWood et al.lll993l ). 
Previous studies also showed that the negative values of the limb-darkening coefficients are 
needed (limb-brightening, e.g. PB94, BPH94) due to same reason of obtaining more precise 
modeling of the light curves. In fact, by adopting negative limb-darkening coefficients for the 
secondary components of UU Sge and V477 Lyr, BPH94 and PB94 achieved more reasonable 
light curve solutions for albedo values greater than one. 

We have performed a number of solutions by taking into account the previous approaches 
to the light curve modeling of the similar systems. The mass ratio 'g' obtained from the radial 
velocity analysis, the gravity brightening exponents of the components {gi=1.0 and (72=0.32), 
and the bolometric albedo of the primary component (y4i=l) were held fixed in the simul- 
taneous solutions. Our first attempt to solve the light curves by using the square-root limb- 
darkening coefficients did not yield a rational solution, thence we decided to use linear limb- 
darkening law. The wavelength-specific linear limb-darkening coefficients of the hot compan- 
ion were adopted as 0.10 and 0. 07 for x^CV) and XifR), res pect ively. These values were as- 



sumed using the tables given by 



Dfaz-Cordoves et al. 



t9m and 



Claret et al. 



bolometric limb-darkening coefficients for both components were adopted from 



9951). Linear 



Van Hamme 



( Il993. note that the limb-darkening coefficients given in these tables are calculated for up to an effective 



Kurucz atmosphere models (lKuruczlll993l ) were used for both stars during the solutions. The 
adjusted parameters were the orbital inclination (z), the mean effective temperature of the 
secondary, (T2-not including the effects that may influence the shape of the light curve such 
as reflection or spots), surface potentials of the components {Qi, Q2), bandpass luminosity of 
the primary component {Li), bolometric albedo of the secondary component (^42), and the 
wavelength-specific linear limb-darkenin g coefficients of the secondary component (x2(V), 
a;2(R) : the input values were taken from I Clare tl (119981 )). The main criterion to secure a rea- 



sonable solution was to achieve a physically realistic value for the albedo of the secondary 
component (0.5^A2^1.0). This criterion allowed us to set a lower limit for the temperature 
of the primary component, Ti, and solution trials were started with a temperature of Ti=73 
000 K. So as to find the best solution, trials were made between the temperatures of 73 
000 and 85 000 K. The temperature of the secondary component was first adopted as 3500 
K considering the expected value of such a low mass star. Iterations were carried out until 



10 M. Af§ar and C. Ibanoglu 



Table 2. Parameters obtained from the sirimltarioons lig'lit curve anah'sifs of UU Sge and V477 Lyr. 



Parameter 


UU Sge 


V477 Lyr 


i (°) 




Q7 1 -i_ n 1 Q 


sn 3*=! + 06 






7 407 -1- f) 040 

1 .'rtc7 ( _1_ W . \J^\J 


13 123 + OQQ 


^2 




3 809 + OOQ 


2 738 + 006 


Ti (K) 




78 000 


49 500 


T2 (K) 




61 Sfi + SO 


3874 + 1 20 




+ L2)v 




QQ40 -1- n n09Q 




+ L2)r 


\J-lJ\JttJ _1_ \J.\J\J£l\J 


n QQ1 1 + fl 00S1 


Ai 




1.0 


1.0 


A2 




0.927 ± 0.037 


0.998 ± 0.036 


91 




1.0 


1.0 


92 




0.32 


0.32 


in) 




0.1423 ± 0.0013 


0.0779 ± 0.0007 


{r2) 




0.2277 ± 0.0010 


0.2091 ± 0.0003 


a;i(bol) 




0.616 


0.616 


X2 (bol) 




0.493 


0.384 


XI (V) 




0.100 


0.192 






0.070 


0.162 






0.478 ± 0.078 


0.210 ± 0.051 


2:2 (R) 




-0.148 ± 0.074 


-0.121 ± 0.053 


rss 




0.007 


0.023 



the convergence was reached at the smallest value for the residual sum of squares (rss). 
Although the solution trials have yielded parameters that are capable of modeling the light 
curves almost equal in appearance, we have achieved the best fit around the secondary min- 
imum for a primary temperature of 78 000 K and a secondary temperature of 6136±80 K, 
where the rss value was at a minimum among other temperatures adopted in the solutions. 
Examination of Ti - rss variation yielded an approximate error of 3000 K for the hot com- 
ponent's temperature. The uncertainty in the temperature of the secondary component, 80 
K, is the formal la error from the solution. The corrected error can be estimated as 240 K 
based on the uncertainty of 3000 K in the effective temperature of the primary. The results 
of the best solution are given in Table 2 (adjusted parameters have standard errors provided 
by the W-D code). The hght curves and the computed fits are plotted in Fig. 3. Some of the 
astrophysical parameters for the components were calculated and are provided in Table 3. 

The mass of the secondary component seems considerably low compared to the mass of a 
main-sequence star with the same temperature. If we assume that the mass of the secondary 
component has not changed dramatically by the end of the CE interaction and it is just a 
low-mass main-sequence star, then the temperature seems significantly higher than that 
expected from an M type dwarf. The possible reasons which could explain this contradiction 
are discussed in § 6. 

Since the simultaneous two-colour light curve solution gives us an opportunity to con- 
strain the temperatures of the components, we believe that the temperatures derived in this 
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Figure 3. V- and R-passband light curves of UU Sge with computed fits. 



paper are mo r e acc urate than the temperatures reported in previous studies. For example, 



Walton et al 



(Il993l ) suggested a prima ry comp o nent t emperature of 50 000 K by comparing 
the merged ultraviolet spectrum with iKuruczl ( 1l99ll ) LTE model atmospheres. They also 
obtained the spectrum of the secondary component during the primary eclipse (with some 
contamination from the sdO star) and assuming that the secondary is a main sequence star 
they classified it as G7 with an uncertainty of three spectral subclasses. BPH94 also took 
some spectra of the cool component during the primary minima. They used atlasg computer 



code (Kuruc2 



19791 ) and generated several atmosphere models to find the best fit to the 



spectrum of the secondary. Under the assumption that the secondary is a main sequence 
star they determined a temperature of 6250±250 K. They used this value to reanalyse the 
light curve of UU Sge obtained by PB93 and gave two different results for the primary com- 
ponent temperature as 57 000±8000 and 87 000±13 000 K, in terms of limb-darkening and 
brightening approximation, respectively. 

The limb-darkening coefficients of t he seco n dary component obtained from the solution 



significantly differ from those given by 



Claret 



( 1l998l ). If we compare our values with pub- 



lished ones, we see that the F-passband linear limb-darkening coefficient corresponds to a 
temperature of ~10 000 K. The i?-passband linear limb-darkening coefficient has a negative 
value and has no reference in the published tables. The only explanation we could make for 



12 M. Af§ar and C. Ihanoglu 

such limb-darkening coefficients would be the reflection effect emanated from the overheated 
face of the secondary, which gives rise to the "limb-brightening effect". 

There ar e several studies concerning the limb-darkening co efficients of ill uminated atmo- 



sphere s (e.g. 



Peraiah 



1983 



Claret 



2004 



Alencar&Vaz 



19991 ). According to 



Alencar fc Vaz 



(119991 ). as the illumination flux increases, the temperature distribution in the atmosphere 
of the secondary component becomes more homogeneous, in other words, center-to-limb 
brightness variation or center-to-limb-darkening over the stellar disk is removed. The ob- 
served increase in the intensity of the brightness of the secondary star from its center t o 
its limb may be explained in terms of horizontal energy transport (IBond fc Graued 119871 ). 
On the other hand, the limb-darkening/brightening coefficients for the systems that have 
strongly irradiated atmospheres still need further investigation. 



4.2 V477 Lyr 



We have analysed the V- and i?-passband light curves of V477 Lyr using the same code 
as described in the previous section. A series of simultaneous solution trials for a large 
temperature range of the primary component were performed. In order to obtain the best 
and realistic solution, we took into account the same criterion for the albedo of the secondary 
component, as we did for UU Sge. This criterion enable us to adopt a starting point for 
the temperature of hot component as 45 000 K. The fixed and adjusted parameters were 
the same as described for UU Sge. Similarly, we ffist attempted to apply the square root 
limb-darkening law during the solutions, as suggested by several authors for the stars have 



temperatures higher than 8000 K (e.g. 



Claret 



19981 ). but we could not reach a rational 



solution and decided to apply the linear limb-darkening law. The input linear limb-darkening 
co efficients for the primary and secondary compone nts w ere taken from the tables published 



by 



Dfaz-Cordoves et al. 



fll995h . 



Claret et al. 



Claret 



(119981 ). Solution trials were 



(Il995h and 

made for the temperatures between 45 000 and 80 000 K. The best solution was achieved 
for the temperature of 49 500 K, where the rss value was at a minimum. And, the error 
in the temperature of the hot component was estimated as 4500 K from the examination 
of Ti - rss variation. The observed and computed light curves are given in Fig. 4 and the 
parameters are listed in Table 2. 

The linear limb-darkening coefficients derived from the solution exhibit a noteworthy 
discrepancy from those published ones that mentioned above. The cool component acts like 



-1.2 
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Figure 4. V- and i?-passband light curves of V477 Lyr with computed fits. 



a very hot star due to the reflection effect caused by the subdwarf primary. The heating on 
the hemisphere of the secondary star that faces to hot primary is such vigorous that the 
irradiation affects the lin ab- darkening over the st ellar surface. Considering the grazing nature 
of the secondary eclipse, iBond &: Grauerl ( 119871 ) made an interpretation on the effectiveness 
of the reflection effect on V477 Lyr and suggested that the horizontal energy transport in 
the atmosphere of the secondary component might be the reas on for the bri ghter limb. 

To determine the temperature of the subdwarf companion, iKalen (119831 ) computed the 
Zanstra temperatures using the H and Hell nebular lines and constrained the primary 
component's temperature in between 32 000 and 60 000 K. In a most recent study, PB94 
presented a new Zanstra analysis based on the H/? and Hell nebular lines and estimated 
th e temperatures as 48 OOP and 8 2 000 K, respectively. They also employed the R-technique 



of 



Schonberner fc Drilling) (119841 ) and derived a temperature of 63 000 ±12 000 K for the 



primary star and considering the temperature range found from several methods they de- 
cided to adopt a temperature of 60 000±10 000 K in the light curve analysis. However, the 
simultaneous two-colour photometry gives us an opportunity to derive more accurate and 
approximate values for the temperatures of the binary components, which we have found as 
49 500 K (fixed) for the primary and 3874±120 K for the secondary star. The uncertainty 
of 120 K is the formal la error from the solution. The corrected error can be estimated as 
350 K based on the uncertainty of 4500 K in the effective temperature of the primary. The 
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Table 3. Astrophysical data for UU Sgc and \' 177 Lat. 



Parameter 


UU Sge 


V477 Lyr 


Ml (Mo) 


0.628 ± 0.053 


0.508 ± 0.046 


Ma (Mo) 


0.288 ± 0.031 


0.145 ± 0.021 


i?l (Ro) 


0.349 ± 0.013 


0.172 ± 0.007 


«2 (Ro) 


0.558 ± 0.021 


0.462 ± 0.019 


Ti (K) 


78 000 ± 3000 


49 500 ± 4500 


T2 (K) 


6136 ± 240 


3874 ± 350 


logpi (cgs) 


5.151 ± 0.050 


5.672 ± 0.053 


log 52 (cgs) 


4.404 ± 0.058 


4.269 ± 0.073 



computed astrophysical parameters are given in Table 3. In fact, the mass of the secondary 
component corresponds to an M type star (as in the case of UU Sge) and a lower temper- 
ature is required for a main-sequence star at the same mass. One of the reasons why we 
could not reach such an expected temperature is the limb-brightening phenomenon, which 
takes place on the overheated hemisphere of the cool components. 

5 THE REFLECTION EFFECT AND THE ENHANCED EFFECTIVE 
TEMPERATURES OF THE SECONDARY COMPONENTS 

The hot subdwarf components of UU Sge and V477 Lyr, which are also responsible for the 
excitation of the nebulae, overheat the hemispheres facing to the substellar points of the cool 
components and the reflection effect reveals itself on the light curves' maxima as a sine-like 
distortion with large amphtudes (Figs. 1 and 2). The irradiation on the secondary compo- 
nents significantly increases the temperature around the substellar points. To estimate the 
reflection rates we applied Fourier analysis to the out-of-eclipse variation of the light curves 
and the following truncated Fourier series was used to compute the Fourier coefficients, 

l^Ao + AicosO + A^coslQ BxsinQ B2sin2e. (1) 

In this equation, / is the normalized intensity, 6 is the phase in degrees, Aq is the normaliza- 
tion level for the intensity, Ai and A2 are the contribution to the reflection effect. Bi and 
B2 are the coefficients represent the other perturbations on the light curve. The normalized 
Fourier coefficients of UU Sge and V477 Lyr for both passbands are listed in Table 4. The 
B coefficients have significantly smaller values than other coefficients, which indicate that 
other perturbations have no remarkable influence on the light curves other than reflection 
effect. 

We have computed the amplitudes of the sine-like distortions caused by the reflection 
effect by using the Fourier coefficients of related systems. The amplitudes for UU Sge were 
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Table 4. Normalized Fourier coefficients (1=1 at 6=90°). 



Targets 


passbands 


Ao 






Bi 


B2 


UU Sge 


V 


1.0265 


-0.1907 


0.0312 


0.0062 


-0.0002 






±15 


±22 


±22 


±15 


±17 




R 


1.0251 


-0.2147 


0.0356 


0.0110 


0.0016 






±13 


±19 


±18 


±12 


±14 


V477 Lyr 


V 


1.0477 


-0.2912 


0.0540 


0.0053 


0.0027 






±22 


±33 


±25 


±30 


±29 




R 


1.0630 


-0.3338 


0.0617 


-0.0011 


0.0015 






±14 


±20 


±17 


±19 


±18 



found to be 0.396 mag and 0.446 mag in V- and i?-passband, respectively. The amplitude 

found in F-passband light curve seems in very good agreement with that of derived by 

PB93. In the case of V477 Lyr, the amplitudes in V- and i?-passband were estimated to be 

0.588 mag and 0.665 mag, in fair agreement with that of obtained by PB94. 

We have tried to estimate the maximum temperature values, T2/1, around the substellar 

points of the secondary components by employing the following equation; 
1 



rp4 

-'-2h 



10 



-0.4Am/ D2m4 



{RiTt + RiTi) - RiTt 



(2) 



Where R2 are the radii and Ti, T2 are the effective temperatures of the primary 
and the secondary components, respectively; Am is the amplitude of the brightness vari- 
ation from minimum (6'=0°, mid-primary) to maximum (0=180°, mid-secondary) of the 
sine-like distortion, and is the approximate central angle viewed from the centre of the 
seco ndary, which corr esponds to the irradiated area over the surface of the secondary star 



see 



Exter et al. 



20051 ). We have determined the V- and i?-passband substellar-point tem- 
peratures of the cool components of UU Sge to be around 30 680 K and 31 800 K, and of 
V477 Lyr to be around 16 860 and 17 550 K, r espectively. These e stimates agree well with 
those of computed by the equation described in Exter et al] (2005). 



6 OVERSIZED RADII OF THE SECONDARY COMPONENTS AND 
EVOLUTIONARY STATUS OF THE SYSTEMS 

Radii of the secondary components of UU Sge and V477 Lyr have been found to be larger 
than those expected for main-sequence stars at the same masses. On the other hand, our 
targets are not the only examples for such systems that contain oversized secondaries. In this 
section, we discuss the possible reasons on the oversized nature of the late-type secondary 
components of our targets along with some other selected PNNi and post-CE systems (listed 
in Table 5). Although the mass of the secondary component of V471 Tau is almost three times 
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the masses of other secondaries studied here , we have espe c ially included this well-studied 



system for comparison. In a study made by lO'Brien et al.l (120011 ). two possible reasons for 



the expanded radius of the K-type secondary of V471 Tau were accentuated: (1) because of 
the possible CE phase, which V471 Tau might have recently emerged from, the secondary 
star may be still out of thermal equilibrium or (2) the starspots covering the surface of the 
star prevent the convective energy from being transported, which as a result increase the 
radius of the star. 

The CE evolution and its effect on the low-mass main-sequence secondary components 
are the phenomena that have been poorly understood both theoretically and observationally. 
If a close binary system, consisting of a WD primary and a red dwarf (RD) secondary, have 
once undergone a CE phase, the mass and the radius of the secondary component might 
have been changed due to mass transfer during the formation of the CE. Although there 
remain questions on whether the mass of the secondary component alters dramatically or 
does not change in a substantial manner, there are observational evidences which show the 
significant changes in the radius of the component, such as the systems given in this paper. 

There are several theoretical studies that have focused upo n the response of the seco ndary 



component to the accretion process during the CE evolution. iPrialnik &: Liviol (Il985l ) made 
some evolutionary calculations for a 0.2 fully convective star which accretes matter at 
several rates and investigated the outcome of the accretion. They found that in the models 
with accretion rates around 10~^-10~^ M© yr~^, the incoming material accumulates faster 
than the contraction rate, and consequently the outer radius increases. In their models the 
final value of the accumulated matter reaches to 2.5x10"^ Mq, which basically indicates 
that the mass of the fully convective star remains almost unchanged. They also concluded 
that the results should be valid for all fully convective stars (^0.3 Mq) as well as for the 
stars with deep convective envelopes (0.3 M0<M<O.5 M0). In a similar study carried out by 



Sarna fc Ziolkowskil (119881 ). the effect of rapid mass accretion onto low-mass main-sequence 



stars (0.3 Mq-O.S Mq) was investigated. They, too, came to a similar conclusion and in- 
dicated that the stars with deep convective envelopes exp eriencing rapid mass accr etion 



significantly increase their radii. Another study published by lHjellming fc TaamI (jl99ll ) also 
probed the response of a 1.25 main-sequence star to the conditions within the CE. At a 
range of relatively high accretion rates (10^^ to 10^^ M© yr^^), they found that the radius of 
the secondary component can expand ten times larger than its main-sequence radius upon 
mass accretion of less than 0.1 Mq. After filling its Roche lobe, the secondary loses back 
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almost all the mass it has accreted and only a small amount of mass (<0.01 Mq) remains 
accumulated after emergence from the envelope. As a result, the surface layers of the sec- 
ondary depart from thermal equilibrium and enter in a thermally relaxation stage with a 
radius larger than its initial main-sequence radius. On the other hand, the results of some 



Sandquist et al 



1998f ) reveal that the 



three-dimensional hydrodynamical simulations (e.g. 
evolution of the CE phase takes place in short dynamic time scales and the secondary star 
which moves rapidly (supersonic) in the CE leaves this stage almost without mass accretion. 
However, the cases we will discuss in the following paragraphs involve the possible small 
amount of mass accretion by the cool companion. 

The second reason which has been considered as a process to increase the radii of the 
low-mass secondary components is the effectiveness of the magnetic activity. There are sev- 
eral studies on low-mass main-sequen ce stars, which compare the physical properties of 



magnetically active and inactive stars. 



MuUan &: MacDonaldl ( 120011 ) showed that, among 



the stars at the same masse s, active stars are co oler and have larger radii than inactive 



stars. In a most recent study 



Morales et al. 



(120081 ) concluded that all active low-mass stars 
(including low-mass eclipsing binaries) are cooler than inactive stars of similar luminosity 
and hence h ave larger radii. 



Recently, 



De Marco et al. 



(120081 ) have discussed the importance of the irradiation effect 
in binary PNNi and suggested that the heating which takes effect on the 'day side' of the 
secondary stars might cause a puffing up and expand the surface layers of the secondaries. 

By taking into consideration two of the reasons summarized above, we have carried out 
some calculations to study the evolutionary status of some selected post-CE systems and 
interpreted which reason may be responsible for the cool components to have oversized radii: 
incomplete thermal relaxation or presence of magnetically active regions. 



KV Vel : The binary nucleus of the PN DS 1. The sdO star was discovered by 



Drilling 



(Il988h and doub. 



Hilditch et al. 



e-lined spectroscopic variations with a period of 8.571 hour were found 



Drilling] (119851). L ater, the renewed results for the system parameters were published by 



(119961). 



BE UMa : It i s an eclipsing bin ary nucleus. Although the light varia 
was first noted by iKurochkinI ( 11964] ) and an eclipse was fir st observed by 



jility of the 



Ando et al. 



the nebula surrounds the binary nucleus was discovered by 



Liebert et al. 



parameters listed in Table 5 and 6 have been adopted from ]Ferguson et al.] (]l999]l . 



1995!). T he system 



system 

km, 



HS 1136+6646 : Discovered as a spectroscopic binary by 



Heber et al. 



(11996] ) and more 
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recently ISing et al.l (120041 ) have indicated that the double-lined spectroscopic system is a 
newly formed post-CE binary. 

PG 1026+002 : The system, which consists of a WD priraary a n d a RD secondary , 
was discovered during the Palomar Green Survey (iGreen et al.lll986l ). ISaffer et al.l (1l993l ) 
reported that the system is a single-lined spectroscopic binary and the secondary compo- 
nent's radius is compatible with it s mass. However, by using the i?-passband light curve 
of the system iBruch fc Diaa (119991 ) found the radius of the secondary as being larger than 
expected from the theoretical models. 

RR Cae : An eclipsing WD-M dwarf binar y system with a period of 7.3 hour. The 



eclipses were first observed by iKrzeminskil ( Il984j ) . iBruchl ( 



Si 



19991) ana 



y sed th e radial velocity 



Bruch fc Diaa (119981 ) and found that 



measurements together with the photometric data of 
the radius of the M dwarf is larger a s compared wi t h a m ain-sequence star at the same 



mass. However, in a more recent study. 



Maxted et al. 



(120071 ) presented ne w phot o metry and 



BruchI (Il999h . the 



double-hned spectroscopy of the system and concluded that, contrary to 
radius of the secondary is consistent with its mass. By taking into account the discrepancy 
between these two results we have performed two separate calculations. 

RE 1016-053 : Post-CE binary consists of a WD pri mary and a RD seco ndar y. Two differ- 



Wood et al. 



(119991 ) and 



Vennes et al 



ent va lues for the secondary star's radius were given by 

(119991 ) depending on two different spectral analysis and the radius of the secondary compo- 
nent was found as slightly oversized. Both radii have been considered here and the results 
are listed in Table 5. 

F^7i Tau : Is a well-known post-CE binary system (WD primary -|- RD secondary) 



and a member of the Hyades star cluster. T he system has a 



observed since the discovery of its eclipses by 



period of 12.5 hour and b eing 



Nelson &: Young 



(Il970h . libanoglu et al.l (ll994J ) 



showed that the RD secondary is an active st ar. Possible reasons 



the secondary component were discussed by 



O'Brien et al. 



br the expanded radius of 



(I2OOII ) and they attributed the 



expansion to the magnetic activity which originates from the starspots partially covering 
the surface of the star. 

We have used the models (age of 1 Gyr, [M/H]=0.0) provided by 



Chabrier &: Baraffe 



(1l997. hereafter CB97I ) to find the "expected" radii, temperatures and luminosities of the 



secondary components. The models have been fitted by sixth-degree polynomials and these 
equations have been used to calculate the ''expected" radii {R2,teo) given in Table 5 (along 
with ''observed" ones). The mass-radius relationship of CB97 is plotted in Figure 5 and 
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Figure 5. Mass— radius relationshi ps for the low mass main— seq uence stars a dopted from different studies. The fits are taken 
from: solid-line: CB97; dotted-line: ICaillault fc Patterso"nl l ll990l) : dashed-line: iLacvl ((iQT^)- Radii of the secondary components 
are plotted as; UU Sge: filled d iamond; V47 7 Lyr: open diamond; KV Vel: cross; BE UM a: asterisk; HS 1136-f6646: open square; 
PG 1 026+002: plus; RR Cae llBruchlll999l'): filled circl e: RR Cae llMaxted et al.ll2007t) : open circle; RE 1016-053 llWood et al.l 
Il999h : filled triangle; RE 1016-053 llVennes et aLlll999l) : open triangle. 
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Figure 6. Mass— effective temperature relationship for the low mass main— sequence stars (CB97). The symbols are the same 
as in Figure 5. 



comparison is made with the models of iLacyi ( 119771 ) and 



Caillault fc Patterson 



fll99nh . The 



ratio of the "ofeserf ec?" and "expected" radii are also given in Table 5. As it is seen from these 
ratios, most of the secondary components, including UU Sge and V477 Lyr, have relatively 
larger radii (2 or higher) than the radii of their main-sequence equivalents. 
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Figure 7. Mass-luminosity relationship for the low mass main-sequence stars. The symbols are the same as in Figure 5. 

The locations of the secondary companions are illustrated in the mass-temperature (M- 
Te) and mass-luminosity (M-L) planes in Figs. 6 and 7. Since the effective temperatures 
could not be determined for the secondary stars of HS 1136+6646 and RE 1016-053, the 
suggested spectral types taken from the referred papers are given in Table 5. In order to 
plot these stars in the M-T p and M-L planes we adopted the effective temperatures from 



the tables given by 



Leggett 



(Il992l ). The temperatures of 3800 K for K7 classification of HS 
1136+6646, 3300 K and 3100 K for M3 and M5 classifications of RE 1016-053 were used in 
our calculations. 

In order to have an idea on the thermal relaxation time-scale (tth) of the envelope 
(accumulated matter) which was formed duri ng the mass trans f er fro m the primary to the 



secondary, we have used tth equation given by 
tth ^ GMM,^,AR/R^L. 



Prialnik fc Liviol fll985h : 



(3) 



Where M and R is the secondary's mass and ^^expected" radius, M^m is the mass of the 
envelope due to mass accretion, AR is the width of the envelope and L is the bolometric 
luminosity of the secondary star. 

Following the theoretical studies mentioned in earlier paragraphs we have assumed three 
Menv values which may have been accumulated by the secondary star: 0.05, 0.01 and 0.005 
Mq. Then we have calculated the tth for each assumption. The results are given in Table 5 



as tth,0.05 5 ^th,0.01 a^Ild tth,0. 



005- 
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Table 5. The observed and theoretical properties of the secondary components of the selected post-CE systems. The description 
of the properties are made in this section. 



Object 


Ma 




T2 


Ta.too 


fi2,too 


-R2,obs/-R2,tco 


tth,0.05 


tth.O.Ol 


*th, 0.005 


Ref. 






(«0) 


(K) 


(K) 


(«0) 




(yr) 


(yr) 


(yr) 




UU Sge 


0.288 


0.558 


6138 


3420 


0.280 


1.989 


1.64x10** 


3.27x10^ 


1.64x10^^ 


this work 


V477 Lyr 


0.145 


0.463 


3874 


3130 


0.170 


2.727 


9.32x10** 


1.86x10** 


9.32x10'' 


this work 


KV Vel 


0.230 


0.402 


3400 


3340 


0.238 


1.688 


1.62x10** 


3.24x10^" 


1.62x10' 


1 


BE UMA 


0.360 


0.720 


5800 


3480 


0.336 


2.142 


1.28x10** 


2.56x10^ 


1.28x10' 


2 


HS 1136+6646 


0.330 


0.790 


K7 


3460 


0.312 


2.530 


2.01x10** 


4.03 xlO'^ 


2.01x10' 


3 


PG1026+002 


0.220 


0.400 


3170 


3330 


0.231 


1.733 


1.85x10^ 


3.71x10'^ 


1.85x10' 


4,5 


RR Cae 


0.095 


0.189 


2500 


2730 


0.118 


1.598 


1.07x10'' 


2.15x10* 


1.07x10** 


6 


RR Cae 


0.182 


0.210 


3100 


3250 


0.202 


1.042 


1.45x10^ 


2.90x10" 


1.45x10'^ 


7 


RE 1016-053 


0.285 


0.357 


M3 


3420 


0.278 


1.283 


4.76x10^ 


9.53x10'' 


4.76x10" 


8 


RE 1016-053 


0.150 


0.265 


M5 


3150 


0.174 


1.521 


2.62x10* 


5.24x10^^ 


2.62x10' 


9 


V471 Tau 


0.930 


0.960 


5040 


5370 


0.832 


1.153 


5.15x10^ 


1.03x10^ 


5.15x10* 


10 



References : The re f erenc es are given for t he ob serve d pr operties of the s econdary components 
(2) Ferguson et al.l lll999h. (3) [Sing et al.l j2004lj. (4) iBm ch & Di; ~ 
iMaxted et al.l ||2007| {. (S nWood et al.l lll999l) . f9') IVennes et al.l II199I 



, 1 ) IffiTditch et al.l lll996t) . 

Dia^ lll999l'). (5) l^f fcr et alJ ||l993l') . (6) iBruchI lll999^ ■ (7) 
(lOl lO'Brienet al.l (120011) . 



Table 6. The properties of the primary components of the selected post-CE systems. 



Object 


Ml 


Ti 


log 31 


Age 


Ref. 




(Mq) 


(K) 


(c9s) 


(yr) 




UU Sge 


0.628 


78 000 


5.151 


3.0x10* 


1 


V477 Lyr 


0.508 


49 500 


5.672 


3.0x10* 


1 


KV Vel 


0.63 


77 000 


5.85 


3.0x10* 


1 


BE UMA 


0.70 


105 000 


6.5 


3.0x10* 


1 


HS 1136-1-6646 


0.63 


70 000 


7.75 


9.0x10^ 


2 


PG1026-I-002 


0.65 


17 600 


8.1 


1.4x10* 


2 


RR Cae 


0.467 


7005 


7.7 


1.9x10^ 


3 


RR Cae 


0.44 


7540 


7.7 


1.9x10^ 


3 


RE 1016-053 


0.57 


55 800 


7.81 


1.8x10" 


2 


RE 1016-053 


0.59 


55 000 


7.84 


1.8x10" 


2 


V471 Tau 


0.84 


34 500 


8.31 


1.2x10' 


2 



References: The references given here are for the models which have been used t o compute th e "age " of the systems (for the 
other properties given in this table please see the references given in Table 5). fl) ISchonbernei] lll983l) , (2) AB98, (3) BA98. 



Since it is very important to know the critical accretion rates, we have also calculated the 
Eddington luminosities of the secondary stars. We have found the critical accretion rates 
to be of the order of a few 10~^ Mq yr^^. Therefore, it is acceptable to consider a mass 
accretion onto secondary components as rapid as a few 10^^ Mq yr^^, which is an adequate 
ac cretion rate leads to inc rease in the radii of the cool components as previously discussed 
bv lPrialnik fc Liviol JlQSsI ). 

For the binary nuclei, UU Sge, V477 Lyr, KV Vel and BE UMa, the presence of the 



PN warrants us in estimating the ages (the ti me since the emergence 



Schonberner 



1983 



Soker 



rom t he CE) of the 
So, we have 



systems as a few 10^ yr (30 000 - 40 000 yr, 
simply adopted the ages as 3x10^ y r for all the nu clei (Table 6). However, we have used 



the evolutionary tracks computed by 



Blockerl (119951 ) . and showed the possible evolutionary 



states of the sdO componets of UU Sge and V477 Lyr in Fig. 8. Although these models 
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Figure 8. The evolutionary status of UU Sge and V477 Lyr. The models are taken from lBlockeil (|l995l V Open and filled circles 
represent the primary component of UU Sge and V477 Lyr, respectively. 



were computed for PNe with single nuclei, the loci of the sdO components seem in fairly 
well agreement (within the uncertainties given for the masses of sdO components) with the 
models have remnant masses of 0.625 and 0.524 Mq, which evolved from the initial ZAMS 
masses of 3 and 1 M©, respectively. In these models, the remnants have CO cores and the 
nuclear energy is mainly provided by the He shell burning. Besides all these approaches, 
we are aware of the fact that our systems are not PNe with "single" nuclei, and binarity 
may have influenced the evolution of the components through the interaction between the 
components and formation of the CE (more details on the evolution of PNe and CE can be 
found in the referred papers), as much as the different mass-loss rates the systems might 
have encountered. 

For the other systems studied here we have u sed the white dwarf evolut i onary models 



from 



Althaus k Benvenutol fll998. hereafter ABQSh and lBenvenuto k AlthausI (11998. hereafter BAQSh 



to calculate the approximate cooling ages, tcooi- In AB98, the evolutionary models were 
computed for the carbon-oxygen WD masses within the range of 0.5-1.2 and surface 
hydrogen layer masses ranging between 10^^'^ and 10^^ Mh/M, where Mis the mass of the 
modeled WD and Mh is the mass of the H envelope. BA98 developed the models for the 
evolution of helium WDs with masses ranging from 0.15 to 0.5 M0, and having outer H 
envelope masses within the range of 10^* ^Al^i/ 4x10^^. And both WD models were 
improved for the metallicities of Z=0.0 and Z=0.001. 
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For the primary masses higher and lower than 0.5 Mq, we have used the tables given 
by AB98 and BA98, respectively. Various Mh/M and both Z values have been taken into 
account and the models have been examined in Te-logg plane. The most applicable models 
have been selected and interpolated for the related primary masses in order to estimate 
the cooling ages of the WDs. We have found that the models for V471 Tau and RR Cae 
having the surface hydrogen layer mass of 10"*^ Mh/M give the best agreement in Tg-log^f 
plane for the metallicities of Z=0.0 and Z=0.001, respectively. However, other systems seem 
consistent with the models have surface hydrogen layer mass of 10"^ Mh/M and metalhcity 
of Z=0.001. Resultant tcooi values are listed in Table 6. 

Since our goal is to search for the more plausible reason on the oversized nature of the 
secondaries, we have compared tth with the ages of the systems as summarized below: 

UU Sge, V477 Lyr, KV Ve/and BE UMa : The thermal relaxation time-scales calculated 
for each of the envelope masses are relatively longer than the ages of the systems, which 
indicate that the surface of the secondary components have not yet reached the thermal 
equilibrium. The discrepancy between the "ofoserf ed" and the "expecteci" radii of the secon- 
daries can be seen in Fig. 5. The difference of the effective temperatures from the theoretical 
ones is also notable (Fig. 6). 

HS 1136+6646 : The comparison of tth and tcooi of the system shows that surface layers 
of the system is still in a thermally relaxation stage which causes the expansion in the radius 
of the secondary. Also the adopted effective temperature which corresponds to a K7 star is 
located above the theoretical fit and appears to be in agreement with this approach. 

PG 1026+002 : When the computed tth values are compared with tcooi of the system, 
depending on the envelope mass, the secondary component seems to be either at the end 
of its thermally relaxation stage or just passed it. However, the low effective temperature 
of the secondary component may be taken as support for the idea that the reason for the 
expanded radius is the presence of the possible magnetically active regions on the secondary 
star. 

RR Cae : For the parameters given by iBruchl (119991 ). the thermal relaxation time-scales 
for the secondary component of RR Cae are shorter than the cooling age of the system. 
The inconsistency both in radius and effective temperature with the theoretical expectations 
indicates the presence of the possible stellar activity. On the other hand, the posibility of the 
higher amount of matter (at least 0.1 Mq) being accreted by the secon dary leaves a question 
mark that the cool component might be still out of thermal equilibrium. 



Maxted et al. 



mm 
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found the radius consistent with the theoretical expectations as seen in Fig. 5. But there 
seems to be a disagreement between effective temperature and its expected counterpart, 
which also supports the assumption of the presence of active regi ons. Be c ause the mass of 
the secondary component is twice as much of the mass found by iBruchI (119991 ) , tth values 
are relatively well below the cooling age of the system. 

RE 1016-053 : There is an interesting result f or the secondary s tar of RE 1016-053. The 



Wood et al 



1999 and 



Vennes et al 



19991 ) 



radius of the cool component (taken from both 
has larger values than expected according to the mass-radius relationship of CB97 and the 
comparison of the thermal relaxation time-scales with the cooling age of the system shows 
that the secondary star still may have thermal instability on its surface layers. But the 
adopted effective temperatures seem to be a little lower than expected, which could tell us 
the possibility of having active regions on the surface of the secondary if we would have 
lower tth values than tcooi- However, our calculations indicate a secondary star which is far 
from being thermally relaxed that makes us to expect higher effective temperatures. Since 
the effec tive temperatures hav e been deterrained from the spectral types estimated by the 



authors flWood et al 



1999 and 



Vennes et al 



19991 ) and they suggest these spectral types re- 



lying on the masses of the secondary star, we believe this is due to probable underestimation 
of either spectral types or the masses. 

¥471 Tau : Our findings on this well-known post-CE system are listed in Table 5 and 
6, but since the mass of the secondary component is higher than the rest of the secondaries 
studied here, it is not plotted in Figs. 5, 6 and 7. It is obvious that the secondary component 
passed its thermally relaxation stage relatively long ago. Nonetheless, its radius seems a 
little larger and the effective tempe rature seems a l ittle l ower. These results are in fairly well 
agreement with those presented by O'Brien et al. ( 200ll ). 



7 THE EFFICIENCY PARAMETER : "acE" 

The CE evolution is one of the most crucial part in the evolution of binary systems. The 
presence of close binaries, for example, consisting of a WD primary and a RD secondary, in- 
dicates that the wide precursor system must have undergone a significant amount of angular 
momentum loss which causes the orbital shrinkage. 

At some stage of the evolution, more massive primary component evolves to AGB (or 
RGB, depending on the initial configuration of the system) and fills its Roche lobe which 
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leads to mass transfer towards the secondary component. If the mass transfer is unstable the 
secondary star can not accommodate the accreted matter at its surface and both components 
are engulfed by the primary's envelope which is later on being called as CE. Since the CE 
does not corotate with the orbiting binary core, the frictional forces arise due to velocity 
difference. This frictional interaction produces gravitational drag that causes components 
to spiral-in together and reduces the binary separation. The change in orbital separation 
results in the change in orbital energy and a fraction of the difference in orbital energy 
is used to eject the envelope from the system by depositing it int o CE. This fraction is 
designated as the '^effjce ngj^ parameter and denoted by "acE" (e.g. iTutukov fc Yungelson 



1979 



Iben &: Livio 



19931 . and references therein) 



bind 



«CE = Xf^' 

where Aii^bind is the difference between gravitational and thermal energy of the ejected 
envelope at the beginning of Roche lobe overflow (RLOF) and AE'orb is the orbital energy 
difference between the beginning and at the end of the spiralling-in process of the binary 
system. 

One of the main goals of this paper is to make a contribution to the studies which 



search for the limits to constrain the value of "qce" parameter. We have mainly 
the appr o ximat ions given by, e.g., [Webbing (120081 ) . 



Han et al. 



fll995h . lDewi fc TaurisI (12000f ). 



bllowed 



De Kooll (Il990l ) in order to estimate the ^^efficency" parameters for UU Sge and V477 Lyr. 



We assume that the secondary components are main-sequence stars and have not been 
affected by the CE interaction. For the estimation of the binding energy of the envelope 
(A_Ebind) we have used the following equation 



bind 



(5) 



where Mi is the mass of AGB star (also assumed as the initial mass of the primary com- 
ponent), Menv = Ml - Mcore, is the cuvelope mass of AGB star, Mcore is the current mass 
of the primary, Rj^ is the Roche lobe radius of the primary (at the onset of mass transfer) 
and A is a dimensionless parameter which depends on the detailed structure of the enve- 
lope and the stellar density distribution. According to Webbink (2007), A is a function of 
me=Menv/Mi=l-(Mcore/^i), and is well-approximated by 

A"^ ^ 3.000 - 3.816me + 1.041me^ + 0.067me^ + 0.136me^ (6) 
The total change in orbital energy is 
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= —2^, 2^0 ' 

where M2 is the mass of secondary component, Oq is the initial separation of the binary 
system, and Of is the final orbital separation. 

In order to calculate these energies one needs to estimate the initial parameters of the 
binary systems. As a first co njecture f or the initial masses of the primaries we have made use 



of the evolutionary tracks of iBlockeii (119951 ). However, by taking into consideration that the 



RLOF may occur at diff erent stages of th e binar y evolutio n we have also used the initial/final 



mass relations given by llben fc Tutukovi (119861 ) (also see llben fc Lividll993l ). If the compo- 



nents experience a case-C RLOF (primary fills its Roche lobe during the He-shell burning) 
then the mass of the remnant is related to initial mass of the primary by Mco~0.5+0.16(Mi 
- 2.3) if 2.3<Mi<4.2 (the masses are in solar units) and by Mco~0.8+0.0875(Mi - 4.2) if 
4.2<Mi<7.6. If the binary system experiences a late case-B RLOF (primary fills its Roche 
lobe shortly before the He ignition. He core is non-degenerate at this point), a degener- 
ate CO or a "hybrid (C+O+He)" remnant develops with an initial/final mass relation of 
Mco~0.167Mi - 0.085 if 2.3<Mi<5, and Mco~0.4+0.07Mi for the range of 5<Mi<10. 

Considering the different stages of RLOF mentioned above, let us first investigate the 
case-C RLOF in which the primary components of UU Sge and V477 Lyr are CO cores. 
Then the initial masses of the primary stars can be estimated as ?»3.1 Mq and «i2.35 Mq, 
respectively. If we assume that the systems have experienced a case-B RLOF and become 
degenerate CO or hybrid cores, then the initial masses of the primary components must 
have been around 4.2 Mq and 3.5 Mq, respectively. By taking these masses as limits for 
the progenitors of the primaries we have performed a series of calculations t o estimate the 



parameters of the systems (following a somewhat similar work to that of 



O'Brien et al. 



2OOII ). 

We used the theoretical AGB core mass-luminosity approximation given by Iben & 
Tutukov (1984, 1985, 1993); 

Lagb = 6 X 10\Mco/ Mq - 0.5), (8) 

where Lagb is the maximum luminosity for a post- AGB star during the plateau phase of its 
evolution and also the luminosity at the time of RLOF. Applying this approximation, Lagb 
values for the current primary masses of UU Sge and V477 Lyr can be calculated as 7680 
and 480 Lq, respectively. 
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Table 7. Possible initial configurations of the systems and the efficiency parameter, ctcE- 



Object 






q 


ao 


Pi 






(Mq) 


{«©) 




(«©) 


(days) 




UU Sge 


3.0 


297 


fO.42 


511 


738 


0.48 




3.f 


294 


10.76 


503 


710 


0.51 




4.2 


266 


14.58 


438 


501 


0.99 


V477 Lyr 


f.O 


60 


6.90 


110 


125 


0.71 




2.35 


45 


16.21 


73 


47 


4.86 




3.5 


40 


24.14 


62 


30 


12.82 


V477 Lyr* 


2.2 


2f0 


15.17 


345 


485 


0.83 



* The sdO component has been assumed as a degenerate "He" dwarf. 



The AGB radii of the p rimaries have been estimated using the AGB mass-luminosity- 
radius relation proposed by Hurley et al.l ( 200ol ): 



Racb = 1.125Mr°-33(L°-^ 



AGB 



0.383r 



0.76 N 
AGB J) 



(9) 



where -Ragb is the AGB radius that the progenitor might have reached at the stage of RLOF. 
The -Ragb values calculated for the different initial masses are also given in Table 7. 

Now Ave can estimate the initial separations, ao, of the systems with the following formula 

jEgdetonlllQSsI ) 

0.49g2/3 

~ 0.6g2/3 + in(l + gi/3)' ^ ^ 

where ri^=Ri^/aQ is the relative radius of the Roche lobe and q is the mass ratio being equal 
to M1/M2. Adopting -Rl=-Ragb the initial separations have been computed for different mass 
ratios (see Table 7). Equations (5) and (6) have been used for the estimation of the binding 
and orbital energies, then the acE parameters have been determined as listed in Table 7. Pi, 
the initial orbital period of the system has been also computed using the Kepler's third law 
and given in Table 7. 

The somewhat problematic point here is that the mass of the sdO component of V477 Lyr 
is of the order of 0.5 Mq which is a critical and controversial value that lies in the transition 
region between the masses of CO and He remnants. This situation gives us a chance to 
consider the sdO component as a degenerate "He" dwarf which once have a progenitor mass 
less than 2.3 M©. If we assume that the primary component filled its Roche lobe right after 
it has de veloped a degene r ate "H e" core, then we can employ the core-mass-radius relation 
given by llben fc Tutukovl (119861 ). 

R = (11) 

Here, R is the radius of primary by the time of RLOF and Mne is the current mass of the 
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sdO star. By adopting this radius as -Rl and assuming that the initial mass of the primary 
was around 2.2 M0 we have calculated the initial parameters and values of Q!ce for the 
system as listed in Table 7. 

8 SUMMARY AND DISCUSSION 

In this paper, we report the results of the first simultaneous two-passband light curve analy- 
sis for the binary PNNi; UU Sge and V477 Lyr. We have also reanalysed the published radial 
velocity data (PB93 and PB94) of both systems in order to increase the sensitivity of the 
orbital parameters employed in the solution of the light curves. The linear limb-darkening 
coefficients were adjusted during the solutions of both systems. Not surprisingly, the resul- 
tant coefficients were found to be inconsistent with the theoretically expected ones, which 
was interpreted as the presence of the limb-brightening effect due to strong irradiation of 
the secondary components. In the case of UU Sge, the best solution was achieved at the 
primary temperature of 78 000 K by yielding a secondary component temperature of 6136 
K. The temperature of the primary component falls into the temperature range of 87 000 ± 
13000 K derived by BPH94. Also the temperature of secondary component indicates a good 
agreement with that of 6250 K estimated by BPH94. 

We have determined the temperatures of the hot and cool components of V477 Lyr as 49 
500 K and 3874 K, respectively. The temperature of the secondary star shows a remarkable 
discrepancy in comparison with the temperature of 5300 K derived by PB94. We believe 
that having two-passband simultaneous light curve analysis with a different approach to 
the solutions yielded more realistic results in terms of the temperature of the components 
of both systems. 

Nevertheless, it has to be stressed that the secondary component temperature of V477 
Lyr, in contrary to that of UU Sge, seems much closer to the value of expected from a 
main-sequence star at the same mass. The discrepancy between the observed and theoretical 
temperatures may be explained in terms of the thermally relaxation stage that the secondary 
stars have been shown to be in. The more the mass is accumulated the longer the thermal- 
relaxation time scale is needed. During the computations, our main assumption was that 
the mass of the secondary stars has not been significantly affected by the CE interaction 
and the highest value for the mass of the accumulated matter was adopted as 0.05 Mq. 
However, we have to bring into attention that the mass of the secondary components might 
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have undergone drastic changes due to accretion, which would also result in the increase of 
the luminosity (and the temperature) of the secondary components. 

The temperature of the primary component of V477 Lyr indicates a good agreement 
with the H/3 Zanstra temperature of 48 000 K while showing a disagreement with the Hell 
Zanstra temperature of 82 000 K, those of derived by PB94. 49 500 K als o lies in the r ange of 
H and Hell Zanstra temperatures of 32 000 K - 60 000 K determined bylKaleJ (119831 ) . PB94 
mentions the possible He ll stratification in the nebula, which leads up to the variety in the 
Hell Zanstra temperatures of the sdO component. The mass loss via thermal pulses may 
explain the likely He stratification and enhanced He abundance (PB94) in the nebula. In 
this case, the Zanstra temperature determined from He ll nebular lines could be misleading, 
which reveals that the temperature determined by using the H nebular lines may carry more 
valuable information on the temperature of the sdO component. 

The comparison of the thermal time scales and the post-CE ages of the systems (Ta- 
ble 5) indicates that the secondary components of the binary nuclei of the PNe are still 
out of thermal equilibrium along with two other post-CE systems: HS 1136+6646 and RE 
1016-053. The rest of the systems that have been reported as having oversized secondary 
components were found to be thermally relaxed relatively long ago. The oversized nature 
and low temperatures of these systems were interpreted as the presence of the possible mag- 
netically active regions which cause the increase in the size of a low-mass main sequence 
star while reducing its effective temperature. Detailed spectroscopic studies of these systems 
are needed. 

Finally, we have implemented several calculations to estimate the common envelope effi- 
ciency parameter, acE, of the systems. We have determined various initial system parameters 
using possible evolutionary scenarios. The possible efficiency parameters for UU Sge were 
found in the range of 0.3 < acE < 1, which is close to the range of previously suggested 
by theoretical studies (e.g. Han et al. 1995). Nevertheless, the parameter values that we 
have found for V477 Lyr, under the assumption of having CO or C+O+He remnant as the 
primary component, indicate a deviation from the theoretical findings. An acE > 1 niay be 
explained as the contribution from other sources (e.g. deposition of recombi nation energy 



Iben&Livio 



i nto C E ejection, nuclear burning on the surface of the secondary star, see 
( I1993I ) for more detail) to the CE ejection beside the efficiency of orbital energy. On the 
other hand, it seems impossible to know which mechanism could have been such effective 
on the efficiency parameter. Therefore, we have found the courage of making a new assump- 
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tion about the primary component of V477 Lyr. By taking into account its current mass of 
0.508 ± 0.046 Mq, we assumed that the primary component is a degenera te He dwarf with 



Iben fc Tutukov 



an ini tial mass of < 2.3 M©. By using the mass-radius relation given by 
(Il986l ). we calculated the radius of the primary at the time of RLOF (independently from 
its luminosity, see equation (11)). These approaches yielded an efficiency parameter of 0.83 
which indicates that about 83 percent of the orbital energy were used effectively during the 
ejection of the CE. 
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